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ABSTRACT 

We measure the spin of XTE J 1550-564 in two ways: by modelling the thermal 
continuum spectrum of the accretion disc, and independently by modeling the broad 
red wing of the reflection fluorescence Fe-Ka line. We find that the spin measurements 
conducted independently using both leading methods are in agreement with one an- 
other. For the continuum-fitting analysis, we use a data sample consisting of several 
dozen RXTE spectra, and for the Fe-Ka analysis, we use a pair of ASCA spectra 
from a single epoch. Our spin estimate for the black hole primary using the continuum- 
fitting method is —0.11 < a* < 0.71 (90 per cent confidence), with a most likely spin of 
a* = 0.34. In obtaining this result, we have thoroughly explored the dependence of the 
spin value on a wide range of mo del- dependent systematic errors and observational er- 
rors; our precision is limited by uncertainties in the distance and orbital inclination of 
the system. For the Fe-line method, our estimate of spin is a* = 0. 55^22- Combining 
these results, we conclude that the spin of this black hole is moderate, a* = 0.49l l20> 
which suggests that the jet of this microquasar is powered largely by its accretion disc 
rather than by the spin energy of the black hole. 

Key words: accretion, accretion discs — black hole physics — stars: individual 
XTE J1550-564 — X-rays: binaries. 



1 INTRODUCTION 

During its principal 1998-1999 outburst cycle, the bright X- 
ray transient XTE J1550-564 produced one of the most re- 
markable flare events ever observed for a black hole binary. 
For ~ 1 day, the source intensity rose fourfold relative to 
neighbouring plateau values, reaching 6.8 Crab. The flux in 
the dominant power-law component rose by t he same factor, 
and then just as quickly its inte nsity declined (|Sobczak et al.l 
l200d : iMcClintock et al.l 120091 ). Four d ays later, AU-scale 
snper luminal radio jets were observed 1 Hannikainen et "all 
120091 ). Their separation angle (~ 255 mas) and relative ve- 
locity (~ 65 mas d _1 ) links the birth of these jets to the 
impulsive X-ray flare. The subsequent detection of large- 
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scale radio jets in 2000 led t o the discovery of relativistic 
X-ray jets l|Corbel et al.ll2002l ). All of the available evidence 
strongly indicates that these pc-scale X-ray and radio jets 
were produced by the unique 7- Crab flare event, and we 
adopt this view. 

The microquasar XTE J1550-564 (hereafter J1550) is 
further distinguished by a pair of high-frequency X-ray 
oscillations with a 2 :3 frequency ratio (184 and 276 Hz; 
iRemillard et al.ll2002T ). During its 1998-1999 eruption, J1550 
displayed all of the active accretion states: hard, steep 
power law (SPL), thermal dominant ( TP) and intermediate 
(INT; IRemillard fc McClintockl 120061 ). The X-ray spectral 
and timing properties of this source h ave been comprehen- 
sively studied by many authors (e.g ., Sobczak et al. 2000l: 
Homan et all l200ll : IRemillard et all 120021 : iKubota fc Pond 



20041 : iDunn et all l201ofi ~ as have the properties of its 
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radio counterpart dCorbel et all l200ll ; IXue et al.l 120081 ; 
lHannikainen et alj|2009f h 

Likewise, the optical counterpart of J1550 was the sub- 
j ect o f a comprehensive dynamical study by lOrosz et al.l 
l|2002t ). The measurement by these authors of a large mass 
function immediately established J1550 as a dynamically- 
confirmed black hole binary with a ~ 10 Mq black hole 
primary in a f .55-day orbit with a late G or early K com- 
panion. This dynamical model was rece ntly revisited usin g 
new photometric and spectroscopic data (|Orosz et alfe oiO). 
Our higher-resolution spectra (60 km s^ 1 ) revealed that the 
mass ratio is extreme (Q ~ 30) and yielded a refined value 
of the mass function, f (M) = 7. 65 ± 0.38 Mr. Of central 
importance to the present paper, lOrosz et alj (|2010h report 
accurate values of the three key quantities that are essential 
for determining the spin of the black hole via the contin- 
uum fitting method, namely the distance D — 4.38lg '41 kpc, 
black hole mass M = 9.10±0.61 Mq, and orbital inclination 
angle i = 74? 7 ± 3? 8. 

Currently, the two principal methods for measuring 
black hole spir Q are modeling th e thermal spectrum of the 
accretion disc llZhang et al.lll997f) and modeling the profile 
of the Fe-Ka line dFabian et al.ll 19891 ; lLaorlll99lh . For both 
methods, spin is measured by estimating the inner radius 
of the accretion disc r- m = Rin/M in standard GR units 
(G — c = 1). i?i n is identified with the radius of the inner- 
most stable circular orbit (-Risco) of the gravitational poten- 
tial and is related to spin via a monotonic mapping between 
the dimensionless ISCO radius -Risco /M and the di mension- 
less spin parameter a* (|Shapiro fc Teukolsky||T983h . Strong 
support for linking R[ n to itjsco is provided by decades of 
empirical evidence that rj n is con stant in disc-dominate d 
states of black hole binaries (e.g., iTanaka fc Lewinl Il995h . 
as shown most compelling ly in our recent stu dy of the per- 
sistent source LMC X-3 jSteiner et alj 1201(f ). Theoretical 
support for identifying Ri n with -R isco is provided by MHD 
simulations of thin accretion discs ^Reynolds fc Fabianl2008l; 



Shafee et all I2OO8I ; iPenna " et al] |2010| ; but see iNoble et all 
20091 . 120101 ). In short, the relationship for thin accretion discs 



between n n , -Risco and a* is the foundation of both the 
continuum-fitting and Fe-Ka methods of measuring spin. 

In the continuum-fitting (CF) method, one determines 
^isco , and hence a, , via measurements of X-ray tempera- 
ture and luminosity (i.e., using X-ray flux, distance D and 
inclination angle i) of the disc emission. In order to obtain 
reliable values of a*, it is essential to (1) select X-ray spectra 
that have a strong thermal component and (2) have accu- 
rate estimates of D, M and i, like those given above for 
J1550. In practice, we fit the X-ray spectrum of the black 
hole's accretion disc to o ur version of the Novikov-Thorne 
thin accretion disc model dNovikov fc Thorndll973l ;lL i et aJj 
120051 ; iMcClintock et al. I I2006T) using an advanced treatment 



of sp ectral hardening (|Davis et al.l 120051 ; iDavis fc Hubenvl 
2006). In this way, we have measured the spins of six other 



1 Black hole spin is commonly expressed in terms of the di- 
mensionless quantity a* = a/M = cJ/GM 2 , where M and J 
are respectively the black hole mass and angular momentum 
llShapiro fc T cukolskv 1983j). Its limiting value is a* = +1 (—1) 
for a maximal Kerr hole rotating in a prograde (retrograde) sense 
relative to the accretion disc; a* = corresponds to a non- 
spinning Schwarzschild hole. 



stellar black hole s. We find spins ranging from a* ~ 0. 1 
jGou et al l lioich to a* > 0.98 jMcClintock et all l2006t ); 
four other spin va l ues are relatively high, a* ~ 0. 7 — . 9 
jShafee et al.ll200d : lLiu~et al.ll2008l . f2010l ; lGou et~aill2009l) . 

In the Fe-Ka method, one determines -Risco by mod- 
eling the profile of reflection-fluorescent features in the disc. 
Most prominent is the broad and skewed iron line, whose 
shape is determined b y Doppler effects, light b ending, and 
gravitational redshift IjRevnolds fc Nowakll2003l ). Of central 
importance is the effect of the redshift on the red wing 
of the line. This wing extends to very low energies for a 
rapidly rotating black hole (a, ~ 1) because in this case 
gas can orbit near the event horizon, deep in the potential 
well of the black hole. Relative to the CF method, mea- 
suring the extent of this red wing in order to infer a, is 
hindered by the relative faintness of the signal. However, 
the Fe-Ka method has the virtues that it is independent 
of M and D, while the blue wing of the line even allows 
an estimate of i. What makes the Fe-Ka method enor- 
mously important is that it is the primary approach to mea- 
suring the spins of supermassive bl ack holes in AGN. The 
spins of several stellar black holes l|Reis et al.l 120091 . l20ld : 
iMiller et ai1l2009l;lBlum et al.ll2009l) a n d supermassive black 
holes jBrenneman fc Reynolds! l2006fc [Schmoll et~al] 120091 : 
iMiniutti et al.ll2009l ; iFabian et al.ll200£ll : IZoghbi et al.ll201C ) 
have been reported using the Fe line method with values 
ranging from a» ~ to a, > 0.98. 

Knowledge of black hole spin has broad importance 
to astrophysics: For example, spin is central to most of 
the many theories of r elativistic jets observed fo r both mi- 
croquasars and AGN (|Blandford fc Znaiek|[l977t l. and it is 
comparably im portant to collapsar models of long GRBs 
jWooslevI Il993l ) and models of black hole formation and 
black hole binary evolution l|Lee et al.l I2002T ) . Hierarchical 
models for the growth of supermassive black holes require 
knowledge of the spin distributions of the mer ging partners 
(|Volonteri et al.ll2005l ; llBerti fc VolonterilbOO&t ). and the ob- 
served properti es of AGN may be strongly conditioned by 
black hole spin dMcNamara et al.ll2009l ; iGarofalo et~alll201oi ; 
iTchekhovskov et alJl201Ch . Spin measurements are likewise 
important to gravitational-wave ast ronomy in predicting th e 
waveforms of merging black holes (jCampanelli et al]|2006l ). 
Knowledge of black hole spin is becoming important to fun- 
damental physics as well, and enlivening questions are being 
asked: e.g., Is the No-Hair T heorem valid and can it be tested 
|johannsen fc Psaltis|[20lbl)? Do we live in a st ring axiverse 
filled with light axions ifArvanitaki et al]|2010h ? 

There have been two prior estimates of the spin of 
J1550. The first of these, a, « 0—0.1, was obtain ed using the 
CF m ethod and a sample of ten RXTE spectra (|Davis et al.l 
2006). This result was based on the old dynamical model 
with approximate val ues of M, i, and D (e.g., D was uncer- 
tain by ~ 45 per cent: IOrosz e t al. 2002). We improve upon 
the work of Davis et al. by using our new dynamical model 
(e.g., with its fourfold better determination of D) and a « 6- 
times larger sample of RXTE spectra, and by our detailed 
treatment of observational and model-dependent uncertain- 
ties. A second 'preliminary' spin estimate, a, ~ 0.76, was ob- 
tained using a spectral model that self-consistently merged 
the disc-continuum and Fe-Ka components, in which the 
spin result for J1550 was predom inately determined by the 
Fe--Ka model ([Miller et al.l |2009| V We will treat the data 
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considered in their study, and additional data; however, un- 
like this exploratory and preliminary study of eight black- 
hole systems (including J1550), we fixedly focus on J1550 
and so are able to refine and improve upon their work. 

In this Paper, we present the spin of J 1550 on two 
fronts. After introducing the data sets (§ [2}, we begin by 
first applying the CF technique (§§ I3l4p . This work is com- 
plemented with a thorough exploration of the statistical and 
systematic uncertainties inherent to the CF method (§ [5]). 
Next, we measure the spin of J1550 using the Fe-Ka tech- 
nique (§ [6}, and we finish with a discussion of the results 
(§[7) and our conclusions (§EJ). 



2 OBSERVATIONS 

The primary data set used in this study is a compendium 
of 347 Rossi X-ray Timing Explorer (RXTE) observations. 
The data include those obtained during the bright discov- 
ery outburst in 1998-1999 on through four additional mi- 
nor outbursts, the last ending in mid-2003. A light curve 
of the flux, showing the spectral evolution of the source, is 
presented in Figure [1] where the spectral state assignments 
have been determined usi ng precisely the model, pro c edure s 
and criteria described in iRemillard fc McClintock] {2006). 
Thermal dominant data, which are of primary importance 
for CF spin measurements, were obtained exclusively during 
the first outburst cycle. 

RXTE spectral data are collected using exclusively 
PCU-2, the best calibrated of the five Proportional Counter 
Array (PCA) detectors. Spectra are individually obtained 
by grouping sequential observations into approximately half- 
day bins with a typical exposure time ~ 3 ks. W e follow the 
procedures described in iMcClintock et al.l (|2006T l: The data 
are background subtracted, a customary 1 per cent system- 
atic error in the data count rates is included, and a dead time 
correction ranging from approximately 1 to 20 per cent is ap- 
plied. Spectra are fitted over the energy range 2.55-45 keV 
using XSPEC version 12.4-1 2.6 (the software used through- 
out this Paper; lArnaud|[l996l ) . Calibration over this range is 
achieved using the latest version of PCARMF (vll.70. 

A linear colli mator correction (ass uming an ideal 1° tri- 
angular response; Ijahoda et al.l 120061 ') has been applied to 
the data set to account for a series of offsets in the PCA 
pointing. Specifically, we normalised the flux upward by 4.4 
per cent (a 2/67 offset) during the 1998-1999 outburst, and 
7.1 per cent (4/28 offset) during the 2000 outburst. After 
April 2001 the correction is just ~ 0.1 per cent. In addition 
to these global corrections, a handful of observations taken 
between 1998 September 7-9 and on 1999 January 6 were 
off-target for unknown reasons by w 0. 2° and required us t o 
make large corrections (for details, sec Stci ner et al.| [2009). 
We estimate that the uncertainty in these flux corrections 
is no more than 1-2 per cent, which has a negligible impact 
on our spectral fitting results. 

In addition to the RXTE data set, we also include a 
25 ks ASCA Gas Imaging Spectrometer (CIS) observation 
taken on UT 1998 September 23 when the source was in 
an intermediate state (see the observation time in Fig. [T] 
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Figure 1. A spectral-state encoded 2—20 keV light curve showing 
all five outburst cycles of J1550. Most disc- dominated data were 
obtained during the primary outburst in 1998-1999. The time of 
the ASCA observation, which is analysed in §[6] is marked at day 
15 by a solid black line. The powerful 7-Crab flare near day 12 
was responsible for ejection of superluminal radio jets. 
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marke d by the black ve r tical line) . Following iMiller et al.l 
(|2009h and IMiller et alJ (|2005r i. we use these ASCA data, 
with twice the resolution of the RXTE PCA data, to exam- 
ine the iron line. However, we do not report a CF analysis on 
these data because the Compton component is too strong. 
We use standard data products for the GIS-2 and GIS-3 
spectra with version 4.0 response matrices. The two spectra 
are fitted jointly over the 1-10 keV energy range, and in 
§ [6] we present our analysis of these data using the Fe-Ka 
method. 



3 CONTINUUM-FITTING ANALYSIS 

We first enumerate our CF data selection requirements and 
define two tiers of data quality. Next, we introduce the first 
and principal of three Comptonised accretion-disc models 
which are applied to the RXTE data set. Then, in the fol- 
lowing section we introduce two alternative models that dif- 
fer principally in their treatment of the Compton reflection 
component. We find very close agreement in the spin esti- 
mate using all three models. 



3.1 RXTE Data Selection 

We first identify two tiers of quality in our data: first- 
class 'gold' spectra are selected from just the strongly disc- 
dominated TD-st ate observations (t he most conducive to 
CF modeling; e.g., IShafee et al1l2006t ). We additionally con- 
sider a broader set of second-tier 'silver' spectra which are 
chosen from either SPL or intermediate states. In order for 
a spectrum to be assigned to either category, it must pass 
the Comptonisation strength, quality, and luminosity filters 
described below. 

In our previou s work on J1550 an d the black-hole can- 
didate H1743-322 ijSteiner et al.ll2009h . we showed that one 
f/pcanrfffflally/s'e a regime of SPL- and intermediate-state spectra 
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to obtain values of n n (or equivalently a*) that are consis- 
tent with those obtained in the (gold) TD state. Based upon 
this work, we select soft thermal spectra with a power-law 
normalisation fsc < 25 per cent, which at fsc ~ 25 per cent 
corresponds to a power-law component roughly an order of 
magnitude stronger than that observed in the TD state. 
Data fulfilling this broader definition (which are not already 
in the TD gold class) are candidates for the silver class. 

To be admitted into either category, we require that 
data meet two additional criteria pertaining to the fit qual- 
ity and disc luminosity. For quality screening, we adopt a 
critical goodness-of-fit, defined as X* /v < 2, and also require 
that the inner-disc radius measurement have a precision bet- 
ter than ri n /ov in = 5. We further require that the luminosity 
of the soft thermal component lie in the range 5-30 per cent 
L-Edd- The upp er luminosity threshold is required by ou r 
thin-disc model (|McClintock et al.1120061 ; IPenna et al.ll2010h . 
The lower threshold eliminates a regime in which the accre- 
tion flow properties may be changing, due to the presence 
of either advective flows I Esin et al. 1997f) or other modes 
of coronal feedback (e.g.. lBeloborodovlll999h . 

3.2 Results I: Continuum Fitting using SMEDGE 

In selecting our data and for determining the 
spin, we employ a variant of the principal model 
that we used in our earlier st udy of J15 5 an d 
H1743-322 (hereafter, Model S; ISteiner et all 120091 ): 

CRABCOR X TBABS X SMEDG E(SIMPLg)KERRBB2) . 

Here we use tbabs (|Wilms et al.l |2000| ') in place of 
phabs to describe the low-energy absorption component. We 
fix the column density to the precise value determined us- 
ing Chandra gratin g data, Ah = 8.0+ ^ j x 10 21 cm -2 (90 
per cent confidence; iMiller et al1l2003l ). The custom multi- 
plicative component CRABCOR simply corrects the response 
of the PCA dete ctor using our standard reference spectru m 
of the Crab (see lSteiner et alj|20ld : pToor fc Sewardlll974l 1. 

Th e key compone nt of this model is ker- 

rbb2 (|McClintock et al.1 l200fj ). a fully relativistic thin 
accretion-disc model, which includes self-irradiatio n of the 
disc ('returning radiation') and limb darkening (|Li et al.l 
2005). The effects of spectral hardening are incorporated 
via a pair of look - up tables for the hard ening factor / 
|Davis et al l 120051 ; iDavis fc Hubenvl 120061) corresponding 
to two representative values of the viscosity parameter: 
a — 0.01 and a = 0.1. Here and throughout, motivated 
by the results of b oth observational data and global 
GRMHD simulations (|Penna et al.l l20ld ; iKing et al.1 120071 . 
and references therein), we adopt a = 0.1 as our fiducial 
value. Following our previous work, we adopt a zero-torque 
inner boundary condition and assume alignment of the 
black hole spin axis with the binary orbital plane; we turn 
on both limb-darkening and returning radiation flags and 
fix the normalisation to one. We fix the input parameters 
M, i and D to their best-fitting values (§ [TJ. The model 
kerrbb2 has just two fit parameters, namely the black 
hole spin parameter a* and the mass accretion rate M, or 
equivalently, n n and the Eddington-scaled b olometric disc 
luminosity Ln/L-Edd <|McClintock et alJlioO^ ). 

We model the power-law com ponent by convolv ing 
the thermal component with SIMPL |Steiner et al.l I2009T ). a 
model that mimics the physics of Compton scattering of 
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Figure 2. The spin, expressed both in terms of R ln and a*, 
versus luminosity. The TD data comprise the gold data set and 
the intermediate and SPL data the silver data set. The mean 
value of R ln is in agreement for the two data sets to within ks 5 
per cent even though the Compton component is much stronger 
for the silver data. 



thermal disc photons by a hot corona. The model SIMPL con- 
verts a fraction fsc of the seed photons into a power law 
with photon index V. We use the standard, upscattering-only 
version. For the reflected component, we assume here that 
the disc elastically backscatters all incident Compton pho- 
tons (generated by SIMPL), apart from a broad iron absorp- 
tion edge feature that is mode lled phenomenologically using 
SMEDGE l|Ebisawa et al.l ll994'). The parameters of SMEDGE 
are the edge energy E E d S e (fitted from 7-9 keV), its optical 
depth r max (unconstrained in the fit), and the width of the 
feature VKEdge (fixed at 7 keV). In the section that follows, 
we consider two models of reflection that are more physically 
motivated. 

Applying our selection criteria to the full spectral model 
yields 35 gold spectra, where most of the winnowing is 
a result of our thin-disc limit on the intrinsic luminosity 
(i.e., prior to scattering) of the accretion-disc component: 
Lo/L^dd < 0.3. We additionally select 25 silver spectra, 13 
of which correspond to SPL-state observations and 12 to 
intermediate-state observations. Our spectral-fitting results 
are summarized in Table Q] (gold spectra correspond to en- 
tries 1-35 and silver to entries 36-60) . 

For all these selected data, in Figure [2] we plot a* versus 
the luminosity of the disc component Lu/L^dd- In a depar- 
ture from our earlier work, in addition to a,, we also plot 
the inner disc radius R ln . The two quantities are equivalent 
in the sense that they are simply related to each other via a 
monotonic analytical formula (§ HJ) . We have chosen to also 
show Ri n because it is the quantity that is more directly 
determined via continuum fitting. 

Figure[2]shows that the gold and silver data sets give re- 
sult that are in good agreement. The net weighted result for 
the combined data set is a, = 0.23±0.07 (r in = 5.22±0.24). 
The gold data give a slightly lower value for the spin, 
a* = 0.20, than do the silver data, a, = 0.27; the corre- 
sponding shifts in the mean value of r- m are respectively +2 
per cent and —3 per cent. 

Although the data are clustered within a few per cent 
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Table 1. Model S Continuum-Fitting Results 
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Note. — 

1. Reported error estimates are symmetric lcr statistical uncertainties. 

2. M, i, and D arc frozen at their fiducial values. 

of a central value of ri n , there is a small ~ 5 per cent in- 
crease with increasing Lo, which is most pronounced for 
Lo /L-Edd > 0.2, a pattern that has been p reviously observed 
for other sources (e.g., GRS 1915+1 05, IrvlcClintock et all 
120061 : LMC X-3, ISteiner et all \201(t ). We tentatively at- 
tribute this effect to a thickening of the disc with luminosity 
and the limitation of our razor-thin disc model. 
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4 CONTINUUM FITTING: TOWARDS A 

SELF-CONSISTENT DISC + REFLECTION 
MODEL 

In the previous section, we used the empirical model SMEDGE 
to crudely account for a prominent spectral feature in the 
reflection component, namely, the broad K-edge of iron. We 
now consider a more physically-motivated treatment of the 
full reflection spectrum, which is generated by that por- 
tion of the power-law flux that strikes the accretion disc 
(I Ross fc Fabianlll993h . We first consider a generalised ver- 
sion of SIMPL that is more appropriate to the problem at 
hand. We then examine two reflection models, ireflect, 
and reflionX, concluding that the former model is better 
for CF fitting, while the latter model is better for fitting the 
profile of the Fe-Ka line (which is considered in §[6]). As we 
describe below, there is presently no unified reflection model 
that is well-suited to both approaches of measuring spin. 



4.1 A Variant of the Power-Law Model SIMPL 

As in § [3] the core of our model consists of kerrbb2 and 
SIMPL. However, we now introduce a modified version of 
SIMPL that is appropriate when including a separate and 
additive reflection component. This model, SIMPL-R, is a 
generalisation of SIMPL that cov ers the two limiting cases 
described by equations 1 & 2 in ISteiner et aD (2009), and 
applies to intermediate cases as well: 



n out (E)dE = (1 - f sc )n in (E)dE 



+ (fsc/x) 



n iu (E )G(E;E )dE 



dE. (1) 



Here, n{ n (E) and n ollt (E) are the seed input and model 
output photon number densities at energy (E). Again, the 
normalisation constant fsc is the fraction of photons di- 
rected into a power law with photon index F, and G(E; Eo) 
is the distribution function of th e output power law (see 
ISteiner et ai1l2009l ; lEbTsawal ll999). The one new parameter 
is x, which determines the fraction of the power-law pho- 
tons that strike the disc. These are the photons which will 
be considered in modeling the reflection component. 

The standard version of SIMPL, which was used in 
the preceding section, assumes that none of the Compton- 
scattered photons strike the disc, or adopting an equivalent 
interpretation, that reflection acts like a perfect mirror with 
no absorption. This corresponds to th e limiting case x = 1 , 
which is described by Equation 1 in ISteiner et al.l (2009). 
In the opposite limit, x = 2, half of the scattered photons 
are redirected downward, illuminating the disc, while failing 
to reach an observer at infinity, As they encounter the disc 
atmosphere, the returning photons are absorbed and repro- 
cessed, thereby generating the re flection component. This 
limit corresponds to Equation 2 in lSteiner et al.l (|2009l ). 

The variant SIMPL-R (Eq.[T} generalises this dichotomy, 
making it possible to treat separately the reprocessed emis- 
sion coming from the illuminated disc via the tunable pa- 
rameter x. This allows one to model a corona quite generally. 
The quantity x— 1 describes the solid angle subtended by the 
disc from the perspective of the corona in units of 2n, which 
we refer to as a covering factor. In this paper, we assume 
that the geometry of the corona is a disc-hugging slab with 



a covering factor of unity (x = 2); thus, half the photons 
escape the system and half strike the disc. As shorthand, 
we will refer to the portion of the Compton component pro- 
duced by SIMPL-R which irradiates the disc as SIMPlC (e.g., 
the second term on the right-hand side of Eq. [T] multiplied 
by the covering factor). 



4.2 Results II: Continuum Fitting using ireflect 
and reflionX 

We first conside r the model ireflect 

l|Magdziarz fc Zdziarskil ll995T l. which computes the re- 
flected spectrum (including scattering and edge absorption, 
but excluding line fluorescence) generated in an ionised disc 
atmosphere that is illuminated by an arbitrary external 
spectrum. We convolve the disc-illuminating component 
SIMPlC with ireflect and isolate the reflected compo- 
nent by setting the parameter reLrefl to -1. (Our model 
implicitly assumes that the observed and illuminating 
power-law spectra are identical.) The ionisation parameter 
£ = L/nR 2 is initially set to f0 4 and allowed to vary freely 
from 10 1 — 10 5 , while the characteristic disc temperature is 
fixed to Tdi sc = 5 x 10 6 K and the metallicity is assumed 
to be solar. Fe-Ka emission is included separately in an 
approximate fashion as an intrinsically narrow Gaussian line 
centred at a restframe energy of 6.5 keV. This composite 
reflection component is then convolved with the relativistic 
smearing kernel kerrconv (|Brenneman fc Reynolds! [2006T ) 
with the radial emissivity index q fixed at the best-fitting 
RXTE value q = 2.5 (see § E}. The complete model, 
which is comprised of an accretion-disc and a power-law 
component, is: crabcorxtbabs(simpl-R®kerrbb2 + 

KERRCONVCE) (IREFLECT® SIMPLC + GAUSS)). 

The primary limitation of this model (referred to here- 
after as Model I) is that alt hough edges are included, the flu- 
orescent line features (e.g., iGarcia fc KallmanlboiOl ). apart 
from Fe-Ka, are missing. Also, the strength of the Fe-Ka 
feature should be tied to the depth of the corresponding edge 
feature, but here that is not possible. Below and in § [5] we 
will demonstrate that these shortcomings of Model I have 
little effect on the CF spin results because for our primary 
gold spectra the reflected component is faint compared to 
the dominant thermal component. However, these issues are 
of critical importance in estimating spin via the Fe-Ka line 
(§©■ 

We now consider a second reflection model, re- 
flionX |Ross fc Fabiani r2005). which we use as a replace- 
ment for (lREFLECT<g>SlMPLC+GAUSS) in Model I given 
above. We will refer to the new composite model as Model R. 
In reflionX, reflection is produced by a power-law spec- 
trum illuminating a cold slab of constant density. The virtue 
of this model is that it properly couples line emission to ab- 
sorption, and it also describes the full Fe-K emission-line 
complex. A major drawback is that it is optimized for mod- 
eling AGN, which have cold discs of lower density. Conse- 
quently, in estimating spin via the Fe-Ka method (§ |6]), 
we use a high-density variant of reflionX, refbhb, which 
inclu des an intrinsic blackbody component (Ro ss fc Fabianl 
l2007t ). Because the blackbody component is hardwired into 
refbhb, it can not presently be used with CF models. We 
therefore use reflionX in concert with kerrbb2 for our 
CF analysis. 
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Figure 3. As in Figure [2] we again plot R ln and a* versus lu- 
minosity, but we now show results for all three of the models 
discussed in the text. The data for Model S, which are repeated 
from Figure [2] show the highest degree of internal consistency. 



In addition to the temperature/density limitations of 
reflionX just mentioned, this model has additional short- 
comings. Of primary importance, it requires that the illumi- 
nating spectrum have a simple power law form. This power 
law is not truncated at low energies and its flux can rival or 
exceed the th ermal flux, thereby leading to unphysical re- 
sults (see e.g.. lSteiner et al1l2009l ). In addition, the strength 
of the reflected component is not linked to the normalisa- 
tion of the illuminating spectrum, and so there is no way 
to ensure that the Compton and reflection components are 
appropriately matched. Nevertheless, we employ Model R 
using reflionX as a second-tier CF model that gives us an 
independent check on the results obtained using Model I. 

In summary, using a variant of SIMPL and considering 
two reflection models, we have progressed toward a model 
featuring a self-consistent treatment of thermal disc emis- 
sion, Compton scattering, and disc reflection. For estimating 
spin via the CF method, we favor ireflect, while for the 
Fe-Ka method we elect to use reflionX and refbhb (§[(>])• 

We now apply Model I (§ I4.2[) to our set of RXTE spec- 
tra, while following the procedures described in § 13.21 In this 
case, we find that only a total of 45 spectra (24 gold and 21 



silver) meet our selection criteria (§ 13. ip . compared to the 
60 selected using Model S. Our fitting results for these 45 
spectra are given in Table [2] 

For the gold spectra, we find excellent agreement be- 
tween the results obtained using Model I, a* = 0.23 ± 0.06, 
and Model S, a* = 0.20 ±0.04. This agreement is illustrated 
in Figure [3l where we also show results for Model R. As is 
apparent, considering only the gold spectra, all three models 
are in excellent agreement - the mean values of r- m are con- 
sistent with one another to within ~ 2 per cent. However, for 
the silver spectra the mean values of ri n are depressed for all 
three models, by ~ 10 per cent for the self-consistent reflec- 
tion models (which track each other closely) and by only ~ 5 
per cent for Model S. Interestingly, the primitive Model S 
performs better than the self-consistent reflection models by 
harmonizing the results obtained from the two data sets and 
delivering the highest degree of internal consistency. Figure[4] 
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Figure 4. Model fits for a gold spectrum (left) and a silver spec- 
trum (right), which correspond respectively to observations made 
on MJD 51121.0 and MJD 51115.3 (see Tabs.[Lj&[2]l. The models 
are differentiated by line colour and the individual components 
by line texture. Note how much weaker the power-law component 
is for the gold spectrum, and how closely all three models track 
the data. 



shows an overlay comparison of the best-fitting results us- 
ing the three models for two representative spectra, one gold 
and the other silver. The total unfolded spectra and their 
components are plotted, as well as the data/model ratio. 
The key result of this section is that using our fiducial val- 
ues of M, i and D ( § HJ) , all three models applied to the gold 
spectra give the same low estimate of spin: 0.15 <a» <0.35. 



5 CONTINUUM FITTING: ERROR ANALYSIS 
AND FINAL SPIN RESULT 

In this section, we broadly consider three sources of observa- 
tional error, both systematic and statistical, which bear on 
our final estimate of the spin. In order of increasing impor- 
tance, these are (1) sensitivity to the details of the spectral 
models employed; (2) X-ray flux calibration uncertainties; 
(3) and the uncertainties in the input parameters M, i and 
D. We then perform a comprehensive analysis that incorpo- 
rates these uncertainties and arrive at our final CF estimate 
of the spin of the black hole. In the following, we present an 
overview; for details, see Appendix [S] 

• Sensitivity to X-ray spectral models. In order to make 
this assessment, we determine the change in r- m when vary- 
ing a single model component or a single one of the parame- 
ter settings. Table [3] gives the mean fractional change in r- ln 
for the gold data sample that arises from changing either a 
model parameter (rows Pl-9) or a model component (Ml-5). 
The 'Change' column lists either the new value adopted for 
the parameter or it describes the change made to the model 
component (where DISC, PL, and REFL refer to the accre- 
tion disc, power law, and reflection components). The third 
and fourth columns respectively list the fractional changes 
in ri n for Model I and Model S. The largest tabulated un- 
certainty arises from the choice of the viscosity parameter: 
Using a — 0.01 instead of the default value (a = 0.1) de- 



8 Steiner et al. 



Table 2. Model I Continuum-Fitting Results 



N 


MJD 


L Edd 








SIMPL-R 








KERRBB2 




IREFLECT 


GAUSS 


xJ/DOF 


State 












r 








fsc 






a* 




M(10 ls g/s) 


£(erg cm/s) 


iV(10 /cm /sj 








1 


51117.4 


0. 


18 


2 


04 ± 





04 


0. 


018 ± 


002 





235 ± 


1)2 1 


3 


37 ± 0.11 


1522 ± 1560 


2.8 ± 2.7 





7/74 


TD 


2 


51119.0 


0. 


If) 


2 


.07 ± 


o 


04 





017 ± 


002 


0. 


230 i 


.02 1 


3 


01 ± 0.09 


1215 db 1345 


2.7 db 2.0 





7/74 


TD 


3 


51121.0 


0. 


14 


2 


06 ± 


o 


06 





006 ± 


.001 


o 


238 ± 


017 


2 


.61 ± 0.07 


1967 ± 2881 


1.4 ± 1.1 





.7/74 


TD 


4 


51124. 7 


0. 


1 1 


2 


07 ± 


o 


Ofi 





012 ± 


002 


o 


240 ± 


.022 


2 


20 ± 0.07 


1179 ± 2224 


2.5 ± 1.5 





.6/74 


TD 


5 


51128.6 


0. 


10 


2 


.09 ± 


o 


05 





020 ± 


003 


o 


236 i 


024 


1 


.89 ± 0.07 


1101 db 1834 


2.9 db 1.5 





.9/74 


TD 


5 


51130.5 


0. 


09 


2 


15 ± 


Q 


.03 





023 ± 


.001 


o 


.221 ± 


018 


1 


.79 ± 0.05 


311.1 ± 114.6 


2.9 ± 0.9 


1 


1/74 


TD 


7 


51132.5 


Q 


08 


2 


.07 ± 


Q 


03 


Q 


018 ± 


.001 


o 


254 ± 


030 


1 


.55 ± 0.07 


1976 db 1868 


2.9 db 1.2 


Q 


8/74 


TD 


8 


51134.5 


Q 


08 


2 


.15 ± 


o 


04 


Q 


013 ± 


001 


o 


.261 ± 


()98 


1 


.43 ± 0.06 


1935 ± 2203 


1.9 ± 0.8 


Q 


8/74 


TD 


9 


51136.9 






2 


08 ± 





y 1 




nru i- n 


008 





276 ± 




1 


.26 ± 0.07 


732.7 ± 5296 


2.0 ± 1.0 


jj 




TD 


10 


51157.6 


n 


12 


2 


.63 ± 







n 




.007 





251 ± 


non 


2 


.33 ± 0.08 


606.7 ± 10500 


0.5 ± 0.8 






TD 


11 


51160.3 






1 


97 ± 









00+0 


.001 





237 ± 


017 


3 


21 ± 0.08 


50600 ± 741400 


0.0 ± 1.5 






TD 


12 


51163.2 






2 


.66 ± 


1 


29 




001 + 


003 





.005 ± 


028 


5 


.62 ± 0.17 


10 ± 3988 


3.4 ± 3.0 




7/74 


TD 


13 


51260.6 




19 


1 


97 ± 





04 






001 





222 ± 


021 


3 


77 ± 0.11 


5823 ± 4795 


4.5 ± 3.0 




4/66 


TD 


14 


51261.8 


0. 


18 


2 


04 ± 


o 


04 





018 ± 


002 


o 


218 ± 




3 


.56 ± 0.12 


1527 ± 1575 


2.4 ± 2.9 





7/66 


TD 


15 


51263.1 






2 


07 ± 





04 






002 





184 ± 




3 


.48 ± 0.15 


1869 ± 1752 


4.7 ± 3.0 






TD 


16 


51264.8 


n 




1 


97 ± 


o 




n 


014 x 


001 


o 


184 ± 




3 


.06 ± 0.12 


5574 ± 3975 


7.2 ± 1.9 




8/fifi 


TD 


17 


51265.6 




14 


2 


10 ± 





08 






004 





219 ± 


020 


2 


.79 ± 0.09 


547 ± 1228 


3.3 ± 1.5 




7/66 


TD 


18 


51266.9 


Q 


13 


2 


08 ± 





23 


Q 


014 ± 


008 





247 ± 


027 


2 


.45 ± 0.12 


672.7 ± 4994 


2.6 ± 2.0 


Q 


5 /66 


TD 


19 


51267.6 




13 


2 


08 ± 


o 






m \ i n 


009 


o 


225 ± 




2 


47 ± 0.13 


697.8 ± 3756 


4.9 ± 2.6 







TD 


20 


51273.6 


n 




1 


.95 ± 







n 




.001 





183 ± 


040 


1 


.97 ± 0.10 


10000 ± 8847 


8.1 ± 1.4 




Q/RR 


TD 


21 


51274.5 






2 


04 ± 











005 





.238 ± 


028 


1 


.68 ± 0.08 


827.3 ± 2862 


3.1 ± 1.6 






TD 


22 


51276.3 




08 


1 


.99 ± 





19 




010 ± 


004 





246 ± 


028 


1 


48 ± 0.07 


516.7 ± 2565 


2.6 ± 0.7 




7/66 


TD 


23 


51278.7 


Q 


06 


1 


92 ± 





06 


Q 


Oil ± 


.001 





341 ± 


031 


1 


11 ± 0.06 


23140 ± 53800 


3.6 ±1.2 


Q 


8 /66 


TD 


24 


51279.6 


— 2. 


06 — 


1 


99 ± 


() 






013 ± 


Oil 


o 


.329 ± 


.022 — 


1 


.07 ± 0.05 


677.9 ± 7036 


2.1 ± 1.1 




_8/66 


TD 


25 


51110.3 


0. 


26 


2 


.45 ± 


() 


02 


— 

0. 


109 ± 


.005 


o. 


233 ± 


.033 


5 


06 ± 0.20 


4879 ± 2066 


16.8 ± 9.5 


- 




8/74 


INT 


26 


51111.6 


0. 


25 


2 


39 ± 


o 


03 





113 ± 


006 


0. 


220 ± 


.040 


4 


.80 ± 0.23 


7147 ± 4072 


24.1 ± 11.0 





.7/74 


INT 


27 


51113.7 


0. 


22 


2 


43 ± 


Q 


03 





108 ± 


007 


o 


234 ± 


.045 


4 


.25 ±0.24 


2331 ± 1531 


16.7 ± 9.4 





.9/74 


INT 


28 


51115.3 


0. 


18 


2 


39 ± 


Q 


07 


0. 


086 ± 


on 


Q 


.317 ± 


.041 


3 


.28 ± 0.26 


875.9 ± 1712 


9.0 ± 6.9 


o 


7/74 


SPL 


29 


51126.6 


0. 


11 


2 


14 ± 


() 


.03 


0. 


025 ± 


001 


o 


.233 ± 


024 


2 


.11 ± 0.08 


10000 ± 6716 


8.8 ± 1.7 


1 


.5/74 


INT 


30 


51269.7 


0. 


12 


2 


.10 ± 





.03 


0. 


022 ± 


001 





225 ± 


025 


2 


.26 ± 0.08 


10000 ± 8200 


9.1 ± 1.9 


1 


4/66 


SPL 


31 


51270.8 


0. 


08 


2 


.17 ± 





.03 


0. 


053 ± 


003 





502 ± 


025 


1 


28 ± 0.06 


10000 ± 6763 


10.2 ± 2.5 


1 


6/66 


SPL 


32 


51271.4 


0. 


08 


2 


.20 ± 





.03 


0. 


059 ± 


003 





469 ± 


044 


1 


.29 ± 0.10 


10000 ± 7370 


9.5 ± 3.0 


1 


2/66 


SPL 


33 


51664.4 


0. 


10 


2 


.40 ± 





04 


0. 


090 ± 


.005 





345 ± 


.029 


2 


78 ± 0.12 


511.2 ± 558.9 


10.2 ± 3.5 





6/66 


SPL 


34 


51664.7 


0. 


15 


2 


.40 ± 





.05 


0. 


097 ± 


007 





407 ± 


.024 


2 


.58 ± 0.13 


568.4 ± 836.3 


10.4 ± 3.8 





9/66 


SPL 


35 


51665.4 


0. 


14 


2 


.37 ± 





13 


0. 


079 ± 


019 





339 ± 


.055 


2 


.49 ± 0.31 


745.6 ± 3070 


7.6 ± 6.3 





6/66 


SPL 


36 


51667.7 


0. 


12 


2 


.29 ± 





03 


0. 


067 ± 


.005 





412 ± 


.029 


2 


.00 ± 0.11 


1533 ± 1302 


8.3 ± 3.9 





7/66 


SPL 


37 


51668.8 


0. 


12 


2 


.32 ± 





04 


0. 


068 ± 


.005 





374 ± 


.034 


2 


.04 ± 0.11 


515.3 ± 602.9 


6.2 ± 2.1 





8/66 


SPL 


38 


51669.2 


0. 


12 


2 


.35 ± 





.02 


0. 


063 ± 


.002 





314 ± 


.023 


2 


27 ± 0.08 


232.8 ± 91.73 


7.6 ± 2.5 





9/66 


SPL 


39 


51670.6 


0. 


10 


2 


.30 ± 





.03 


0. 


079 ± 


.004 





423 ± 


.023 


1 


72 ± 0.07 


455.3 ± 353.1 


8.0 ± 2.3 





7/66 


SPL 


40 


51670.8 


0. 


10 


2 


.25 ± 





.03 


0. 


060 ± 


.004 





439 ± 


.030 


1 


.68 ± 0.09 


1661 ± 1319 


7.2 ± 3.2 





7/66 


INT 


41 


51671.4 


0. 


11 


2 


.29 ± 





13 


0. 


058 ± 


016 





340 ± 


.036 


1 


.96 ± 0.17 


646.3 ± 2684 


7.3 ± 2.5 





9/66 


SPL 


42 


51672.4 


0. 


09 


2 


.25 ± 





10 


0. 


067 ± 


014 





439 ± 


.040 


1 


.53 ± 0.15 


818.8 ± 2652 


7.4 ± 4.4 





7/66 


SPL 


43 


51673.0 


0. 


10 


2 


.33 ± 





07 


0. 


109 ± 


013 





504 ± 


.045 


1 


.49 ± 0.16 


947.3 ± 1858 


10.1 ± 5.7 





6/66 


INT 


44 


51673.4 


0. 


10 


2 


.36 ± 





.09 


0. 


103 ± 


017 





440 ± 


.060 


1 


.61 ± 0.23 


860.7 ± 2522 


10.1 ±6.1 





8/66 


INT 


45 


51674.7 


0. 


08 


2 


22 ± 





.03 


0. 


113 ± 


006 





556 ± 


.046 


1 


.20 ± 0.11 


2605 ± 1802 


11.3 ± 5.6 





6/66 


INT 



Note. — 

1 Reported error estimates are symmetric la statistical uncertainties. 

2 M, i, and D are frozen at their fiducial values. 

creases n n by ^ 6 per cent for Model I and 3 per cent for 
Model S. Each of the other 13 changes considered affect n n 
by < 3 per cent for either model. 
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• Flux calibration. The problem of flux calibration is en- 
demic to X-ray astronomy. Th e Crab spectrum, as deter- 



Table 3. Systematic Changes to the Model 



mined bv lToor fe Sewardl (l974h . is the widely-adopted stan- 
dard that we have consistently used in our work. Uncertain- 
ties in the nor malisation of this spectr um have recently been 
considered by IWeisskopf et~ai1 (|2010h . Using their Figure 1 
as a guide, we adopt a generous ±10 per cent uncertainty in 
our overall flux calibration, which corresponds to a 5 per cent 
uncertainty in r- m . 

• Uncertainties in M, i and D. As in our earlier work 
fe.g.. lLiu et"aill2008l ; [Gou et al.ll2009l , l2010l ). we sample the 
allowed parameter space assuming Gaussian errors (except 
here for D, we use an asymmetric Gaussian). The sampling 
is performed using 42,500 triplets of M, i, and D, which 
are distributed in a uniform grid throughout the parameter 
space. At each point in the grid, the complete RXTE data 
set is analysed with Model S, and the selection criteria given 
in § El are separately applied to the results. Folding all of 
the runs together, a composite distribution based on all of 
the selected spectra is obtained, where we have addition- 
ally weigh ted over the set of possible dynamical models (see 
Table 1 in lOrosz etaHl2010h . 

In conducting this analysis, we have included the robust 
no-eclipse constraint, i < 82°. We have further required that 
during the TD-state plateau phase (days 105-182; Figure[T]) 
the disc luminosity not exceed 85 per cent of LEdd (the 
ac tual Eddington limit fo r disc geometry; see Section 6.1 
in iMcClintock et al.ll2006l ). Lastly, we also require that the 
disc luminosity during the thermal plateau phase be greater 
than 10 per cent of LEdd, or else the full sample of TD data 
would extend downward in luminosity to the unreasonably 
low value of < 0.1 per cent I/Edd- 



In the analysis described above, we have used the de- 
fault value of the viscosity parameter, a = 0.1. Because a 
is the major source of uncertainty considered in Table [3] we 
have repeated the analysis just described using a = 0.01 and 
combined the two distributions, weighting them equally. We 
combine all other errors in Table [3] yielding an ensemble 
value of m 4.2 per cent. Finally, we add in quadrature the 5 
per cent error in the absolute flux calibration and arrive at 
our net error of 6.5 per cent. The effect of this uncertainty on 
our measurement of spin is incorporated by running a box- 
car smoothing kernel (with a 13 per cent full width) over 
the distribution for n n - 

The dominant source of error is the observational un- 
certainties in M, i and D, which in turn are largely due to 
the uncertai nties associated wi th modeling the optical/NIR 
light curves l|Orosz et alJl201Ch . Figure [5] shows the depen- 
dence of 7"i n /spin on these model parameters. Here, using 
the results of the grid analysis described above, we vary one 
of the three parameters, fixing the other two at the their 
best values. The strong correlations between spin and in- 
clination, and between spin and distance, demonstrate the 
degree to which measurement errors in these quantities con- 
tribute to the uncertainty in spin. Together, errors in M, 
i, and D account for Aa» ~ 0.25 (Ari n /r- m ~ 0.2) for the 
90 per cent confidence interval. The contribution due to the 
inclination is sizable, ~ 11 per cent, because its uncertainty 
and its value are large (< 4 degrees and 74.7 degrees, re- 
spectively). The uncertainty in D likewise contributes ~ 11 
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Figure 5. Probability contours for the relationship between 
In /spin and M, i, and D. Each of the first three panels (a-c) 
shows variation for a single parameter; the other two parame- 
ters have been fixed at their best values. The orientations of the 
probability ellipsoids show that spin is positively correlated with 
M and negatively correlated with both i and D. In panel (d), a 
combined probability distribution for the case a = 0.1 is shown 
(arbitrarily scaled) with variation in M, i, and D folded together. 



per cent, while uncertainties in the mass affect the spin only 
at the ~ 7 per cent level. 

After folding together all sources of error, the result- 
ing probability distribution is shown in Figure [6] with R ln 
and a* displayed respectively on the bottom and top axes. 
The green vertical line identifies the most probable spin, 
a» = 0.34, and the yellow lines indicate the 90 per cent 
confidence interval, which extends from -0.11 to 0.71. From 
an inspection of this distribution function, we conclude that 
the black hole is unlikely to be in a retrograde configuration 
(only ~ 11.2 per cent probability). Of greater importance, 
we conclude that the spin is not high. For example, the prob- 
ability that the CF spin exceeds 0.9 is less than 0.4 per cent, 
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K,„ (M) 

Figure 6. Composite probability density for R; n and a* which 
takes into account both systematic and statistical errors includ- 
ing: uncertainties in distance, black hole mass and inclination; 
the spectral model; and the uncertainty in the absolute flux cal- 
ibration. The net probability distribution is a combination of 
the individual distributions for two values of a. The contribu- 
tion from each integrates to 50 per cent probability and is shown 
for a = 0.01 (blue dash-dotted line) and a = 0.1 (red dashed 
line). The 90 per cent confidence limits for the combined distri- 
bution are shown as yellow vertical lines, the lcr limits as vertical 
black lines, and the most probable spin is marked with a green 
line. We conclude that the spin is moderate: —0.11 < a, < 0.71 
(90 per cent confidence). 



a surprising result for a black hole that has produced super- 
luminal jets. 



6 SPIN FROM REFLECTION FEATURES 

In the previous section, we concluded that the spin param- 
eter has a low or intermediate value. This result is based 
on our CF analysis of many RXTE spectra, which were 
obtained primarily in the TD state. In what follows, we 
first analyse intermediate-state spectra of J1550 obtained 
with the ASCA X-ray observatory; we then supplement this 
analysis using a sample of RXTE spectra, also obtained in 
the intermediate state. We fit the ASCA GIS-2 and GIS-3 
data simultaneously, using a floating normalisation constant 
to allow for cross-calibration uncertainties. Our work dif- 
fers from_earlier analyses of these same data by others (e.g., 



Miller et al.ll2005l : iGierliriski fc Pond 120051 : iMiller et all 



120091 ): Our focus is on a detailed analysis of the reflection 
component, rather than on a precise model of the overall 
continuum. We begin by setting all the physical parame- 
ters of the binary to the best-estimate values presented in 
lOrosz et al l (|201QT i. 

In addition to the soft disc and hard power-law com- 
ponents seen in the TD- and intermediate-state spectra of 
black hole binaries, a bro ad emission line at ~ 6.4 keV is 
also often present (see e.g. Miller 2007). This line feature is 
merely the most prominent reflection signature that arises 
as hard emission from the corona irradiates the cooler disc 
ijRoss fc Fa bian 2005). In the vicinity of a black hole, the 
lvoxi-Ka line shape and other refl ection features are dis- 
torted by various relativistic effects (|Fabian et al ]|l989l ; lLaorl 




Energy (keV) 

Figure 7. Data/model ratio for a phcnomcnological continuum 
model consisting of a thermal disc and a Compton component. 
ASCA GIS-2 and GIS-3 spectra are shown in black and red re- 
spectively. The data were fitted jointly in the 1—4 and 7-10 keV 
energy range. The residuals in the 4-7 keV band show the rel- 
ativistic nature (asymmetry and broadness) of the iron-emission 
line profile. The data have been rebinned for plotting purposes. 
The data/model ratio for the full energy range is shown in Fig. [8] 



Il99lf ). The spin parameter can be constrained by modeling 
these features because their shape depends on how far the 
disc extends down into the gravitational potential well (see 
§[l]), the key assumption again being that this extent is set 
by the radius of the ISCO. 



6.1 Phenomenological Models - ASCA 

In order to highlight the relativistic nature of the line pro- 
file in the ASCA spectra, we start by modeling the 1— 
4 keV and 7-10 keV continua with a combination of a 
disc blackbody (des cribed by the XSPEC model DlSKBtjf] of 
iMitsuda et~al]|l984h and the Comptonisation model SIMPL. 
The neutral hydrogen column was i nitiall y fixed at Ns = 
8 x 10 21 cm" 2 as per IMiller et ahl l|200ot ), which resulted 
in a poor fit to the continuum with x 2 /v = 2172.3/1002. 
Allowing the column density to vary resulted in a signifi- 
cant improvement to the fit with \ 2 l v ~ 1367.5/1001 for 
iV H = 5.4 ± 0.1 x 10 21 cm~ 2 . The total neutral hydrogen 
column density in the line of sight to J1550, which was de- 
termined using the Chandra transmission gra ting, is not ex- 
pected to vary (|Miller. Cackett fc Reisl [2009). However, the 
differing values of Nn can be reasonably attributed to dif- 
ferences in the calibrations of the Chandra and ASCA de- 
tectors. Furthermore, allowing A^h to differ between the 
two GIS spectra further improves the fit: A\ 2 = —47.5 
for one less degree of freedom with a difference in TVh of 
< 5 per cent. Figure [7] shows the line spectrum obtained 
by modeling the continuum as described above. The asym- 
metric and broad residual feature in the 4-7 keV band has 
the appearance one expects for fluorescent disc-line emission 
arising near a black hole. 



3 This model characterises the thermal emission using only two 
parameters — the flux normalisation and a colour temperature. 
Here, we use this very approximate model of the continuum (com- 
pare §[3} because of its simplicity in phenomenologically describ- 
ing the thermal continuum. 



The Spin of the Black Hole in XTE J1550-564 11 



We provide a physical description of the Fe line by 
first m odeling the residuals seen in Figure [7] using the LAOR 
model (|Laorlll99lh and fitting for the inner radius n n and 
the power-law index q of the emissivity profile, which is de- 
scribed by a power-law of the form e(r) oc r~ q . The outer 
disc radius is fixed at the maximum allowed value of 400 r g 
(r g = GM/c 2 ), and the disc inclination is const rained to be 
appro ximately la from the adopted value of lOrosz et al.l 
( |2010l ) (i.e. between 71 and 78 degrees). The line energy 
is constrained between 6.4-6.97 keV. The fit achieved by in- 
cluding the LAOR component, shown in Figure [51 results in 
X 2 /u = 1848.0/1501, an improvement of A X 2 = -416 for 5 
fewer degrees of freedom (compared to the best-fitting con- 
tinuum model with no line feature). The best fit parameters 
for this model are detailed in Table [4] (Model 1). 

It can be seen from the ratio plot shown in Figure[8jthat 
this simple, heuristic model, although mostly adequate, does 
not provide a detailed description of all the features present 
in the 6-8 keV range. Adding a narrow Gaussian line at 
~ 6.7 keV only marginally improves the fi10 (Ax 2 = -14.4 
for 2 fewer degrees of freedom), with evidence for additional 
residuals, which are possibly associated with Fe- if -shell ab- 
sorpt i on edges in par tially ionised material (|Ross fe Fabianl 
1 19931 : iRoss et all 1 19961 ). Such features are usually present 
at ~ 7.1 keV in TD-state spectr a of black hole binaries 
|Done et al.ll 19921 : iReis et alj|2008l ). In order to properly ac- 
count for the panorama of features associated with the re- 
processing of radiation in the accretion disc, we now consider 
complete reflection models. 
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Figure 8. (top:) The ASCA spectra. Below are plots of the ratio 
of the data to a phenomenological continuum model consisting of 
thermal-disc and Compton components plus a LAOR line; in the 
bottommost panel, a narrow Gaussian line has been added to the 
model. 



6.2 Reflection Analysis - ASCA 

We replace the la or component with reflionX 
l|Ross fe Fabianl 120051 ). which describes the spectrum 
reflected from an optically-thick and cold atmosphere 
of constant density that is illuminated by a power-law 
spectrum (§H]). The parameters of the model are the iron 
abundance (set to Solar), photon index of the illuminat- 
ing power law, ionisation parameter, and normalisation. 
The gravitational and Doppler effects are accounted for 
using the fully relativistic co nvolution model kerrconv 
ijBrenneman fe Reynolds 2006), which includes black hole 
spin as a fit parameter. The power-law indexes of reflionX 
and the Compton component (simpl-R) are tied, and, as 
before, we constrain the inclination to lie between 71 and 78 
degrees and include a narrow Gaussian line at w 6.7 keV. 
The model results in a good and improved fit to the data 
with x 2 /v = 1752.3/1499 (Model 2 in Table QJ; however it 
still does not fully account for the reflection features, with 
residuals present at m 7 keV (top panel of Fig. [9]). 

Although the scattered fraction for this spectrum is 
high, fsc > 50 per cent, and the CF method is not ap- 
plicable, we nevertheless investigated the effect of switching 
the continuum model from diskbb to kerrbb (Model 3), 
with the mass, distance and inclination frozen at their nom- 
inal values. This change produced insignificant differences in 
the fit parameters (Table [4]). For both Models 2 and 3, we 



4 This feature was previously associ ated llTom sick et al.l |200lf) 
with emission from the Galactic ridge iValinia fe Marshall 199^) . 
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Figure 9. Data/model ratio for (top:) the reflionX model to- 
gether with a separate thermal emission and Compton compo- 
nent, (bottom:) Self-consistent thermal emission and reflection 
(refbhb) together with a power-law component. 

find that the spin parameter is moderate (< 0.75). Mean- 
while, the disc ionisation has pegged at its maximum value 
(£ = 10 4 erg cm -2 s _1 ) indicating that the surface layer of 
the accretion disc is highly ionised, with iron possibly being 
fully ionised. In such circumstances, the Fe absorption edge 
can be parti cularly strong an d is often found to be highly 
smeared (see lRoss et al.|[l996l and references therein). 

In order to incorporate the effects of thermal ionisation 
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Table 4. ASCA 1-10 keV spectral fit parameters with a variety of reflection-based models. 



Model 1 Model 2 Model 3 Model 4 Model 5 



Parameters SIMPL0DISKBB SIMPL-R®DISKBB SIMPL-R(g)KERRBB POWERLAW 

+LAOR ±KERRCONV(g)REFLIONX +KERRCONV(gSREFLIONX +KERRCONV0REFBHB 



N H (xl0 22 cm" 2 ) 


0.576+g'gg 

— U.UUz 


0.650+2!™ 

— U.UU.J 


— U.UU f 


0.663+g-gg 

— U.UUO 


0.653 n nn7 


r 


2.40 ±0.01 


9 qoq+0.006 
z.ozy_ Q 01Q 


9 Q9n+0 ° 03 
z.ozu_ 002 


2.24 ±0.01 


9 99+O.O3 
z,z -0.02 


fsc (N hard ) a 


0.6 ±0.3 


0.64 ± 0.04 


0.616 ±0.002 


2.4 ± 0.1 


9 q+0.2 
z,a -0.1 


kT (keV) 


n crio+0.009 
u,ol -0.006 


566+0- 001 
u.ooo_ 013 




0.540 ±0.001 


542 +0.002 


Ndiskbb ((n/iofcpc' aos *) 


52Q0±«g 


421 1 + 86 








M (xlO 18 g s" 1 ) 






0.668+ ,™ 3 






q 


2.0 ±0.2 


1.88 ±0.01 


1 85+ ' 2 


9 Qo+0.04 


9 c+0.2 


i (degrees) 


71-78 


71-78 


71-78 


77 ±1 


82_3 


E L aor (keV) 


6.40+0.01 










r-m Og) 


8-2±l;l 










Nlaor (XlO" 3 ) 


7.1 ±0.1 










£ ( erg cm s _1 ) 




10000- 320 


10000_9oo 






Nreflionx (xlO -6 ) 




1.32 ±0.06 


1 oso+0.004 
1.2»8_ 07Q 






H den (xl0 22 H cm" 3 ) 








1.00-0.02 


1.00-0.02 


^illum 1 * 6b 








99+0- 03 


95+0- 08 
u ' zo -0.07 


ATre/bhi, (Xl0~ 2 ) 








6 2 +a2 


5.96+°J 


spin (a*) 




< 0.75 


0.45(< 0.75) 


0.6(> 0.38) 


U.OO_ 22 


x 2 /^ 


1848.0/1501 


1752.3/1499 


1759.5/1499 


1700.9/1498 


1698.6/1498 



Notes: All errors are quoted at the 90 per cent confidence level for one parameter of interest A\ 2 = 2.71). Model 1, which is purely 
phcnomcnological, uses the familiar LAOR line and allows a comparison with previous work. Models 2 and 3 use different disc 
components; however both of them employ the same full reflection model (reflionX), while treating the Compton component using 
SIMPL-R (§ 14.11 1. The core of Model 4 is REFBHB which is likewise a full reflection model, with the added virtue that it self-consistently 
models the thermal component as well. In Models 1 to 4 the inclination was constrained to be between 71 and 78 degrees. In Model 5 
the inclination is allowed to range from 60 to 82 degrees. A constraint on the inclination was achieved only for Models 4 and 5. 
a The POWERLAW normalisation is in photons cm -2 s — 1 for Models 4 and 5. For Models 1-3, the normalisation is given by the 
dimensionless parameter /sc (see §|3j. 



expected for a hot accretion dis c, we replace reflionX with 
the model refbhb developed bv lRoss fc Fabianl (|2007l ) . This 
reflection model accounts for both thermal X-ray emission 
and the reflection features. The effects of Compton broad- 
ening in the disc are fully included, subject to the one as- 
sumption of a constant-density atmosphere. The parame- 
ters of the model are the number density of hydrogen in 
the illuminated surface layer, ifden, the temperature of the 
blackbody heating the surface layers, the power-law photon 
index, and the ratio of the total flux illuminating the disc 
to the total blackbody flux emitted by the disc. Again, we 
tie the power-law index of refbhb to that of the Compton 
component - now modelled as a standard power law - and 
convolve the spectrum with kerrconv in order to include 
relativistic broadening. The model results in an excellent fit 
to the data with % 2 l v = 1700.9/1498 (Model 4, see bottom 
panel in Fig. however the hydrogen surface density is 
pegged at the maximum value of the model. The ionisation 
state of the disc is inversely proportional to the value of the 
hydrogen density and thus the pegged value implies that the 
fit is requiring a higher amount of emission in the form of 
discrete features as opposed to the near featureless reflected 
continuum arising from a highly ionised disc-surface. A sim- 
ilar result would be produced by increasing the iron abun- 
dance. Unfortunately the model in its current format does 
not allow for a change in elemental abundances. In order to 
investigate the effect that Hden has on the spin parameter 



fixed it at 1 x 10 21 H cm 3 using Model 4 (i.e. an or- 
der of magnitude less than the value presented in Table [3} 
and refitted the data. This constraint on Hd En resulted in an 
adequate fit with v 2 /v — 1731.6/1499 and a spin value of 
0.60 ± 0.05. We note here that powerlaw has been used to 
model the Compton component. We have explored replacing 
powerlaw with SIMPL-R, and the fit becomes worse with 
X 2 /v = 1802.5/1498. However, the value of the spin param- 
eter, as well as those of the reflection parameters, remains 
largely unchanged. 

From Models 1 to 4 it is clear that the spin parame- 
ter is consistently below ~ 0.75. However, in the first three 
cases inclination is not constrained. For this reason we ex- 
plore a very broad range of i, from 60° to 82°. We note that 
above this limit, the disc would be super-Eddington dur- 
ing its steady thermal plateau in Figure [T] The best fit is 
given by Model 5 (Tab. [4| and reaches the upper inclina- 
tion limit, netting a small improvement (A^ 2 = —2.3) over 
Model 4. For all models we see that the emissivity index is 
consistently below the typical value of three associated with 
the canonical 'lamp-post' coronal geometry, and is instead 
more consistent with a slab-like corona. In order to illumi- 
nate any degeneracy between the value of spin and either 
the emissivity index or inclination, we show in Figure [10] 
the 68, 90 and 95 per cent probability contours for these 
parameters plotted versus spin. In both instances there ex- 
ists a small and negative correlation with spin. However it is 
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Figure 10. (left): Emissivity versus spin contour plot for J1550. The 68, 90 and 95 per cent confidence range for two parameters of 
interest are shown in black, red and green, respectively. We have allowed i to take any value between 60 ^ i ^ 82, and find that the spin 
is greater than 0.33 at the 90 per cent level of confidence, (right:) Similar plot for inclination versus spin. We see from the ASCA data 
that a zero spin value is clearly ruled out, as is an inclination lower than 72°. 
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Figure 11. Goodness-of-fit versus spin parameter for J1550. 
From the reflection features present in the ASCA spectra of J1550 
we can rule out a non rotating black hole at over 3<r confidence. 
However we cannot place a comparable strong upper limit on 
this value. The 90 per cent confidence level (A% 2 = 2.71 for one 
parameter of interest) is shown in magenta. The black dotted lines 
indicate confidence intervals. Spin is constrained to 0.33 < a* < 
0.70 at the 90 per cent confidence level. 



also clear that q is well constrained between 2.2 and 2.7 and 
that i > 75° at 90 per cent confidence even while includ- 
ing the uncertainty in spin. When we marginalize over these 
parameters and compute the uncertainty in spin alone (Fig- 
ure [11] for Model 5), the spin parameter obtained from the 
gravitational blurring of reflection features is constrained to 
be in the window 0.33 < a, < 0.70 at 90 per cent confi- 
dence with the best estimate at a* w 0.55. A non-rotating 
Schwarzschild black hole is rejected at greater than 3a. 

Our measured spin using the refbhb model is consis- 
tently lowe r than the prelimin ary value of a* = 0.75 — 0.80 
reported bv lMiller etafl (2009), and we have attained a bet- 



ter fit than they (Ax 2 < -100) for more degrees of freedom. 
The critical difference in our model and spin estimate comes 
from having incorporated the effect Compton-broadening of 
the iron Ka line in the hot layers of the accretion disc. 
With refbhb the disc is intrinsically hot and therefore the 
effect of Compton-broadening is fully accounted for when 
modeling the data. The extra broadening caused by this ef- 
fect acts to lower the degree of gravitational broadening and 
as such requires less extreme spin parameters as compared 
with models where the reflection is as sumed to come from a 
relatively cold surface, e.g, for AGN l|Ross fc Fabianl [20071 1. 
This behavior can indeed be appreciated when we compare 
the spin value obtained from reflionX - a reflection model 
specifically designed for a cold accretion disc - to that of 
refbhb. We see from Table Q] that for Models 2 and 3 the 
spin cover a higher range, with the 90 per cent error going 
as high as 0.75. 



6.3 Spin from reflection features — RXTE 

In order to supplement the ASCA spin measurement above, 
we present an analysis of a sample of ten RXTE spectra 
selected from the composite data set discussed in §§ 13141 to 
have the following properties: very large scattered fraction, 
fsc > 50 per cent, goodness of fit, \ 2 l v < 2, and uniform 
values of luminosity and photon index, Lo /L^dd ~ 0.2T0.05 
and r w 2.5 ± 0.1, respectively. We begin by simultane- 
ously modeling the reflection features present in all the 
RXTE spectra using reflionX convolved with kerrconv 
while using simpl-R®diskbb for the thermal plus Comp- 
ton continuum (Model 2 in § I6.2f) . (N.B. The refbhb com- 
ponent used in Models 4 and 5 was unable to converge to 
an adequate fit for RXTE and provided no spin constraint. 
Therefore, in this section we adopt Model 2.) The spin, in- 
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clination (60 ^ i ^ 82) and emissivity indejjfl are treated 
as global parameters among the ten spectra. As in § [3] the 
neutral hydrogen column density is fixed at 8 x 10 21 cm" 2 . 
The remaining parameters were allowed to vary in individ- 
ual spectra. Figure [12] shows the best-fitting model spectra 
(top panel) together with the data-to-model ratio for each 
spectrum (bottom panel). The fit was marginally improved 
(Ax 2 = -30.5 for 10 de grees of freedom) by including a 
narrow line at «6.7 keV which accounts for the slight cur- 
vature in the residuals at that energy (compare the lower 
two panels in Fig. I12[) . 

We find that the global best fit is sensitive to the up- 
per energy range adopted for the RXTE spectra, which we 
attribute in part to a competition between the lower en- 
ergy reflection features and the high-energy Compton hump. 
Considering upper ranges between 12-45 keV, the best spin 
estimate was found between a, ~ 0.6 — 0.69 giving reduced 
chi-square values from x 2 jv — 0.4 — 0.8, with higher values 
obtained at extended energy ranges. Most importantly, the 
model consistently estimated the 90 per cent upper limit for 
spin at a, — 0.75. For the other global parameters, we treat 
the RXTE results as second-tier, but find results consistent 
with the ASCA values: q w 2.5 and i > 72° (90 per cent). 

We are cautious in interpreting this spin estimate using 
RXTE spectral fits, owing to the coarse (~ 20 per cent) 
energy resolution. However, we expect that RXTE should 
provide robust upper bounds on the degree of relativistic 
broadening (viz., spin), owing to it s vast collecting are a and 
< 1 per cent spectral calibration ( Jah oda et al.l I2006T I . We 
caution towards the significance of the RXTE -derived spin 
parameter and consider the upper limit obtained here as a 
complementary result to that obtained from the ASCA data 
alone, confirming that Fe-Ka spin is not high. 



7 DISCUSSION 

7.1 A Combined Fe-Ka and CF Result 

In the two previous sections, we concluded that both the Fe- 
line and CF methods predict moderate values of spin, which 
are quite consistent: 0.33 < a» < 0.70 (ASCA only) and 
—0.11 < a» < 0.71, respectively (90 per cent confidence). 
The CF spin result predicts a slightly narrower Fe-line fea- 
ture than that found by the Fe-line analysis. Alternatively, 
the Fe-line measurements consistently favor a high inclina- 
tion, and therefore require a lesser distance (D «4 kpc), in 
order for the CF results to match. 

Having obtained two independent measurements of the 
spin, we now combine them by convolving the individual 
spin probability distributions to obtain the joint distribu- 
tion shown in Figure 1131 Our synthesized result is then 
0.29 < o» < 0.62, with a most probable value of a* = 0.49. 
Remarkably, based on a model of bi n ary ev olution and the 
GRB collapsar model, iBrown et al.l (|2007l ) predicted that 
J1550 formed with a, w 0.5. We confirm their prediction. 



5 Because we have selected a homogeneous set of spectra with 
almost identical luminosities, it is likely that the emissivity index 
- an indicator of coronal geometry - is the same for all ten spectra. 
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Figure 12. Best-fitting model (top) for the ten RXTE obser- 
vations. The model consists of a disc and Compton continuum 
together with a relativistically blurred reflection component. The 
spin, inclination, and emissivity index were treated as global pa- 
rameters (see § 16. 31 1. (middle:) Data/model ratio for the above 
model and (bottom:) after the inclusion of a narrow Gaussian 
emission line. 
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Figure 13. Combined Fe-Kct and CF probability density for R ln 
and a* . The net result is again a moderate value in-between the 
two individual estimates (a* = 0.49^q 2q, 90 per cent confidence). 



7.2 Confronting GRMHD Simulations 

Recently, it has become feasible, via general-relativistic mag- 
netohydrodynamic (GRMHD) simulations, to assess the dif- 
ferences between MRI-driven accretion flows and the ide- 
alised cv-disc model, upon which our CF model is based. 
These differences include both a non-zero torque inside the 
ISCO and an al tered angular mome ntum profile for the disc. 
The results of iPenna et alj (|2010j) show, for the geomet- 
rically thin accretion discs we consider, that these effects 
are expected to impact our results by < 7 per cent in r ln 
(roughly Aa* < 0.09 here). In general, they conclude that 
deviations from the Novikov-Thorne model are small for thin 
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discs at low luminosity and grow worse as the luminosity 
rises and the disc thickens. 

Similar MHD simulations have been made to assess 
a principal assumption of the Fe-line method, namely, 
that the line emission fro m within the ISCO is negligible 
^Reynolds fe Fabianll 2008 ) . Including the effect of contribut- 
ing plunging-region emission results in intrinsically broader 
line profiles and hence will lower the estimate for spin. For 
the disc thickness and spin values in question, simulations 
predict that this effect could possibly shift n n by ~ 12 
per c ent (value taken from Fig. 5 in I Reynolds fc Fabianl 
1200a ) , thereby decreasing the most probable Fe-Ka estimate 
of spin from a, m 0.55 to m 0.4, in very close agreement with 
the best CF value of o» « 0.34. 

We infer that both spin estimates are subject to mod- 
erate corrections. For the Fe-Ka reflection method, we can 
specifically conclude that the corrected value of spin is lower 
than the measured value, strengthening our conclusion that 
the spin of J1550 is moderate. 



7.3 The Question of Alignment 

The spin of an accreting black hole in a binary is expected 
to align with the orbital angular momentum vecto r of the 
system within « 10 r - 10 s years (iMaccaronei 120021 ). A re- 
cent population synthesis study, which makes conservative 
assumptions concerning the torques acting to align a black 
hole, predicts that most black holes will be aligned to bet- 
ter than 10° . In § [6] we constrained the inclination of the 
inner, reflecting portion of the accretion disc (Fig. 1 10p . This 
allows us to check on the relative alignment of the black 
hole spin axis (which is aligned with this inner-disc region; 
lLodato fc Pringla I2006T I and the orbital vector. In our ex- 
ploration of the Fe-Ka model, for a wide range of orbital 
inclinations (60°- 82°), we find a best-fitting inclination for 
the inner disc of « 75 — 82° . This value is consistent with 
the orbital in clination angle giv en by our dynamical model, 
i = 74.7±3.8 (lOrosz et alj|20ld), whi ch validates the CF as- 
sumption of alignment ( Li et al]|2005h . while simultaneously 
providing support for the dynamical model. 



7.4 Implications of a Low-Spin Microquasar 

The low spins of J1550 and the microquasar A0620-00 
(a, « 0.1; iGou et alJlioioh challenge the long-standing and 
widely-held belief that there is a stron g connection between 
black hole spin and relativistic jets (Blandfor d fc Znaiekl 
1 19771 ; hereafter BZ). We note that this belief is also chal- 
lenged by a statistical study that found no e vidence for a 
link between black hole spin and jet power dFender et al.l 
I2010T ). In any case, if jets are powered by black hole spin, 
then theory predicts that jet power will increase d rarnat- 
ically with increasing a* i|Tchekhovskov et al.|[201tih . with 
the jet receiving more power fro m the accretion disc than 
from the black hole for a* > 0.4 (|McKinnevll2005l ). 

Given the low spins of J1550 and A0620-00, it would 
appear that their episodic jets are driven largely by the ac- 
cretion disc. One well-known candidate mechanism is the 
ce ntrifugally driven outflow o f matter from a disc described 
by iBlandfoird fc Pavnelll982l (hereafter BP). A useful com- 
parison of the operational regimes of BP and BZ is given by 



iGarofalo et"aH (|2010h . They show that BP is always viable, 
but that BZ is a more likely mechanism for the most rapidly 
rotating sources, such as the extreme- Kerr BH microquasa r 
GRS 1915+105 (|McClintock et al.ll2006l ; iBlum et alj|2009h . 

The relativistic, two-sided jet of J1550 was launched 
during the remarkable 7 Crab flare (see § [TJ . We now show 
that during this daylong event the luminosity of the accre- 
tion disc was close to, or perhaps at, its Eddington limit, 
which suggests that radiation pressure was a key agent in 
collimating or feeding a magnetically-accelerated jet via a 
disc win d (for a discussion of the interplay between jet and 
wind, see lNeilsen fc Leell2009l and lMiller et al.ll2008h . In Fig- 
ure fTJ] we plot the intrinsic accretion-disc luminosity during 
the flare state versus the luminosity during the thermal- 
dominant plateau state (Days 105-182; see Fig. [T]). Each 
data point represents a single triplet of values of M, i and 
D from among the 42,500 triplets considered in our Odyssey 
cluster analysis (Appendix A. 4), and that data point was 
derived by analysing the complete data set for J1550. The 
spin for each point (averaged over the gold and silver data) 
is indicated by a colour, which is encoded in the bar at the 
top of the plot. Th e point correspond ing to the dynamical 
model adopted by lOrosz et al.l (|2010h . Model F, which is 
our fiducial model, is labelled and marked by a red cross. 
The five less probable models considered by Orosz et al. are 
marked by black crosses. We conclude that Model F, by far 
the most probable model (see Appendix A. 5), is very near 
the Eddington limit for disc accretion. 

Luminous discs are geometrically thick and widely be- 
lieved to be effective at driving jets. In the case of the J1550 
flare, the disc is not only thick, it is also near its Eddington 
limit, so that it will provide substantial radiation pressure 
and possibly even shed material via a radiation-driven out- 
flow, thereby promoting a jet-ejection event. In any case, as 
a bottom line, the low spins of J 1550 and A0620-00 indicate 
that spin is not the sole driver behind all powerful episodic 
jets. 



8 CONCLUSION 

For the first time in a single work, we have determined high- 
quality estimates of the spin of an accreting black hole using 
the two independent, leading methods. In our CF analy- 
sis, we carefully explored the sensitivity of our results to a 
wide range of model-dependent systematic errors and ob- 
servational errors. We conclude that J1550 is a slowly spin- 
ning black hole with a» ~ 0.34, while ruling out spins larger 
than a* > 0.71 at 90 per cent confidence. Next, we analysed 
the Fe-Ka and reflection signatures in bright, intermediate 
spectral states of J1550. By modeling these broad, skewed 
features, we obtained a slightly higher estimate of the spin, 
o* « 0.55lg'22 (at 90 per cent confidence), while also de- 
riving an estimate of the inclination angle of the inner disc 
that is in close agr eement with the orbital inclination angle 
l|Orosz et al.l I201CJ ) . Combining the two spin estimates, we 
conclude that J1550, like the microquasar A0620-00, is a 
slowly spinning black hole. 

The low spins of both J1550 and A0620-00 indicate 
that, for at least some microquasars, BZ-type mechanisms 
are not primary in driving powerful episodic jets, and that 
other mechanisms (perhaps BP) are at play. The near 
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Figure 14. The intrinsic (i.e., seed) luminosity of the disc com- 
ponent during the 7 Crab flare versus the luminosity during the 
thermal plateau phase. In order to avoid saturating the plot, we 
show only half the data points, which were selected at random. 
The vertical black lines mark the lower and upper luminosity 
thresholds, and the horizontal dashed line corresponds to the Ed- 
dington limit of an accretion disc (§0. Note that Model F is very 
near this limit. 

Eddington-limited 7 Crab flare observed for J1550 suggests 
that radiation-pressure support from a thermal disc is one 
possible way that low-spin black holes are aided in driving 
large-scale relativistic jets. 



We allow the column density (P2-3) to vary over a 
broad range, Ah = 6— 10 x 10 21 cm -2 , which corresponds to 
> 8cr relative to the precise value determined using Chandra 
grating data (see § 13. 2f) . We consider this extreme range 
because of the discrepant results for Ah obtained using 
ASCA data (see § O, which we attribute to an error in 
the calibration of the ASCA detectors at low energies. As 
shown in Table [3] our liberal estimate of the uncertainty in 
Ah affects our determination of n n by < 3 per cent. 

We next explore the parameters of the ireflect model. 
We test smaller covering factors of 1/2 and 1/10 by linking 
the covering factor in ireflect to —(a; — 1) from the SIMPL- 
R model, where the leading minus sign acts as a switch in 
the model to isolate reflection from the direct (illuminating) 
component. Thus, we consider two cases: x = 1.5 (P4) and 
x — 1.1 (P5). We next try fitting for the covering factor, 
allowing it to vary between and 1, while also fitting for 
the emissivity, which we constrain to lie in the range 2 < 
q < 5 (P6). As shown in Table [3] the effect of this exercise 
on ri n is small w 1 per cent. Smaller still is the effect of 
varying the disc temperature. Decreasing Tdi sc (P7) by a 
factor of 5 relative to its assumed value (§ 14.21) . we find 
that the ionisation parameter increases slightly, but that the 
effect on n n is negligible (< 0.1 per cent). 

Lastly, we adjust the width of the SMEDGE component 
WEdge to first half (P8) and then twice (P9) its nominal 
value of 7 keV, which impacts ri n by < 0.5 per cent. 

In summary, as we found earlier in our study of LMC 
X-3 (|Steiner et alJl201Ch . a is the parameter (aside from M, 
i, and D) that introduces the largest uncertainty in deter- 
mining spin via the CF method. 
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APPENDIX A: CONTINUUM-FITTING: 
ASSESSING THE SYSTEMATIC 
UNCERTAINTIES 

Al Model Parameters 

We consider the effect of the principal parameters listed in 
Table[3]on our final determination of the spin (Figure[6]) for 
Model S (§EJ) and Model I (§[3]). Here and below, we consider 
only the gold data. As in §§ 13.214.21 we fix M, i, and D at 
their fiducial values (§[!}. As illustrated in Figure [6l the 
effect of decreasing the viscosity parameter (PI in Table [3]) 
is to decrease n n (by ~ 3 — 6 per cent, depending on the 
model). 



A2 Model Components 

We begin by substituting bhspec |Davis fc Hubenvl I2OO6D 
for the thermal disc componen t in place of kerrbb2 (see 
§ 14.21 in iMcClintock et aT]|2006l for a discussion of these rel- 
ativistic disc models and their relationship). The virtue of 
bhspec relative to kerrbb2 is that it directly computes the 
effects of spectral hardening; its drawback is that it does not 
include returning radiation, which heats the disc. Employ- 
ing bhspec and following the procedures described in the 
preceding section, we find that r- m is increased by ~ 1 — 3 
per cent (Ml in Table©. 

Next, we explore the possibility that the power-law com- 
ponent is cut off exponentially at high energy (e.g., thermal 
Comptonisation), while allowing the cut-off energy to vary 
over the range kT e = 25 — 200 keV. (We did not correct 
fsc in order to achieve photon conservation because this 
correction is negligible for the gold spectra.) We find that 
the effect of a possible cutoff is small, changing n n by < 1.5 
per cent (M2). In addition, we generated the power law us- 
ing the double-sided version of SIMPL and SIMPL-R (in place 
of the upscattering-only version). The effect on n n is < 0.3 
per cent (M3). 

Lastly, we examine the effects of substituting one of 
the reflection components for another. We find that both 
ireflect and reflionX give somewhat smaller values of 
9"i n than SMEDGE, but the effect is small, < 3 per cent (M4, 
M5). 
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A3 Flux 

As described in § [5] we include a liberal ~ 10 per cent al- 
lowance for the uncertainty in the absolute-flux calibration. 
Because the luminosity of the thermal component at a given 
colour temperature scales proportionally to r ; 2 n , a 10 per cent 
adjustment to the flux normalisation introduces a 5 per cent 
uncertainty in r[ n . 

A4 Black Hole Mass, Inclination and Distance 

While analysing the X-ray spectral data, we have used the 
best estimates for M, i and D fj [1]) taken from to Model 
F of Table 1 in lOrosz et~a l. ( 2010). To explore the sensitiv- 
ity of spin to uncertainties in these measured values, we use 
the Odyssey computing cluster at Harvard University and 
fit the data set at each point in a 3-D grid of 42,500 points 
distributed uniformly over mass, inclination, and distance. 
The grid spans the ranges M = 5 — 17.5 Mq, i = 36 — 85° 
and D — 3 — 7 k pc, respectively. We adopt the 3 kpc distance 
bound following lHannikainen et alj l|2009t) ; the 7 kpc bound 
is a relat ivistic limit based on the proper motion of the X- 
rav jets ilCorbel et al]|2002l) : D ^ c/*/JJQm- < 7 kpc (e.g., 
iMirabel fe Rodriguezll 19991 ). At each grid point, we compute 
a tab le of the spectral hardening factor (e.g., see iGou et al.l 
l20ld ) and fit all of the available TD, SPL, and INT spectra. 
We use Model S (§ I3.2[) because it is computationally effi- 
cient and has the best performance of all three CF models 
considered. We perform the analysis for both values of disc 
viscosity: a — 0.01 and a = 0.1. 

We then apply our data selection criteria, obtaining a 
sample of gold and silver spectra (typically 50-100) at each 
of the 42,500 grid points. From this we derive a spin prob- 
ability distribution unique to each point. Before summing 
over the grid, we impose the following grid point selection 
constraints: First, the grid point's i nclination must be below 
the e clipsing limit, i < 82° (see e.g-. lNaravan fe McClintockl 
l2005h . Secondly, as discussed in § [5] we require that the 
intrinsic disc luminosity during the TD-state plateau phase 
(days 105-181; Figure[fj| fall in the range 0.10 < L D /L Edd < 
0.85. We combine the distributions for all satisfactory grid 
points, weighting each according to its location in the grid 
(with high weights occurring at probable values of M, i, and 
D). 

A5 Rolling Together the Uncertainties 

We combine the systematic uncertainties discussed above in 
two stages. Referring to TableO in the first stage we combine 
in quadrature the individual values in the Model S column 
for rows P8-9, Ml, M2-3, M4-5 with half the value for P2-3 
(half because the range of variation considered for Nn is so 
extreme). For each of the parameters Nh and W^Edge, we use 
the larger of the deviations given in the table. The resultant 
error of 4.2 per cent is combined with the 5 per cent error in 
rin from flux uncertainty to give a net error of 6.5 per cent. 
This combined uncertainty sets the half-length for a boxcar 
smoothing kernel that we apply to the full spin distribution. 

The resulting distribution is shown in Figure IA1I Be- 
cause we have so far considered just dynamical Model F, 
the distribution of r- m is narrower, —0.14 < a, < 0.57 (90 
per cent confidence level) than our final distribution shown 
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Figur e Al . Sim ilar to Figure [6] but using just Model F in 
lOrosz et al.l ll2010h . The green, black, and gold vertical lines in- 
dicate the most likely value for R ln (a*), and 1<t and 90 per cent 
confidence interval limits, respectively. 



in Figure [6] although the most probable value of spin is 
unchanged, a* = 0.34. We now go on to the second stage 
in combining sources of error and consider an ensemble of 
possible dynamical models. 

The case of J1550 is unusual in that there are sev- 
eral candidate models which produce reasonable fits to 
the dynami c al dat a, which are summarized in Table 1 in 
lOrosz et al.l l|2010h ). Above, we considered only Model F, 
the most probable model. We now incorporate the possibil- 
ity that one of the five alternative models (Models C-E and 
Models G & H) are correct. Models A and B do not con- 
strain the dynamical model satisfactorily, and do not allow 
one to obtain a useful distance estimate, and so they are 
disregarded here. 

As was done above for Model F, a spin (n n ) proba- 
bility distribution is obtained for each candidate dynamical 
model. We use for each model, including Model F, the total 
X 2 (summed for the velocity data and the light curve data, 
both optical and infrared) to determine its corrected Akaike 
Infor mation Criterion (AICc; lAkaike|[l973 : Irlurvich fc Tsail 
Il989h . which is closely related to the log-likelihood of each 
model. Using these values, AlC-weights are assigned to 
each model (i): Waic, i = Ex pf-1/2 (AICc, - m/{AICc})] 
IIBurnham fc Andersodl2002h . 

Our fiducial dynamical model is by far the most likely, 
carrying ~ 84 per cent of the total weight. A weighted sum is 
computed using the AIC- weights to obtain a composite spin 
distribution. This is broadened using the boxcar smooth- 
ing kernel described above (13 per cent width) to produce 
the final distribution as shown in Figure [U Thus, this final 
result incorporates uncertainties in the choice of the dynam- 
ical model; dynamical model uncertainties; the X-ray spec- 
tral model and model parameter settings as summarized in 
Table 0- and a 10 per cent uncertainty in the X-ray flux 
calibration. 
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